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ABSTRACT 

A previously inconspicuous knot in the M87 jet has undergone a dramatic outburst and now exceeds 
the nucleus in optical and X-ray luminosity. Monitoring of M87 with the Hubble Space Telescope 
and Chandra X-ray Observatory during 2002-2003, has found month-timescale optical variability in 
both the nucleus and HST-l, a knot in the jet 0.82" from the nucleus. We discuss the behavior of 
the variability timescales as well as spectral energy distribution of both components. In the nucleus, 
we see nearly energy-independent variability behavior. Knot HST-l, however, displays weak energy 
dependence in both X-ray and optical bands, but with nearly comparable rise/decay timescales at 220 
nm and 0.5 keV. The flaring region of HST-l appears stationary over eight months of monitoring. We 
consider various emission models to explain the variability of both components. The flares we see are 
similar to those seen in blazars, albeit on longer timescales, and so could, if viewed at smaller angles, 
explain the extreme variability properties of those objects. 

Subject headings: galaxies: individual (M87) — galaxies: active — galaxies: jets — BL Lacertae 
objects: general — magnetic fields — radiation mechanisms: nonthermal 
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1. INTRODUCTION 

M87 is among the nearest galaxies with a bright 
radio/optical/X-ray jet. Its proximity (distance 16 Mpc, 
Tonry 1991, giving a scale of 78 pc/arcsec) allows fea- 
tures to be studied with unparalleled spatial resolution. 

The dynamic nature of the M87 jet was first recognized 
in the X-rays by Harris et al. (1997, 1998), who found 
variability on 1 year timescales. In the optical, Biretta 
et al. (1999) found superluminal motion throughout the 
inner 10" of the jet, with speeds up to 6c, as well as 
gradual changes in fiux. These observations did not ex- 
plore timescales ^ 6 months, however. More recently, 
the jet's optical spectrum was found to harden in knots, 
consistent with local particle acceleration (Perlman et 
al. 2001a). And Chandra observations found spectral 
indices a^; > 1 {Si, oc broadly consistent with an 

extrapolation of the radio-optical synchrotron emission 
(Marshall et al. 2002, Wilson & Yang 2002). X-ray 
synchrotron emission implies radiative lifetimes ~ 1-10 
years, and thus requires in situ acceleration. 

All this evidence suggested the possibility of shorter- 
timescale variability in the M87 jet. Indeed, recently 
Harris et al. (2003, hereafter paper I) detected X-ray 
flaring on timescales ^ 1 month during 2002 in both 
the nucleus and HST-l, a knot in the jet. Here we dis- 
cuss observations of the M87 jet with both HST and the 
Chandra X-ray Observatory, which find fiares in both the 
nucleus and knot HST-l on timescales of 2 months. A 
second paper (Biretta et al. 2003, hereafter paper III) 
discusses the longer timescale variability of the M87 jet. 
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2. OBSERVATIONS AND DATA REDUCTION 

We make use of Chandra and ifST observations of M87 
during November 2002 - June 2003. During this interval 
we observed M87 roughly every 20-40 days. Table 1 lists 
the _ff5'Tobservations; the Chandra data will be discussed 
in an upcoming paper (Harris et al., in prep.). 

Our Chandra data reduction procedures were detailed 
in Paper I. All HST observations were reduced in IRAF 
and PyRAF with the best available flat fields, biases, 
darks and illumination correction images. Dithered im- 
ages were combined using PYDRIZZLE, which com- 
bines drizzling (Fruchter & Hook 2002) with cosmic ray 
and hot pixel removal, geometric correction and mosaic- 
ing. Galaxy subtraction was done on the F475W and 
F814W data using ELLIPSE, BMODEL and IMCALC. 
Identical galaxy models were used at all epochs. Flux- 
calibrated images were obtained using SYNPHOT. We 
applied a standard Galactic extinction law and assumed 
N{H) = 2.4 X 10^°cm~^ to correct the fluxes for redden- 
ing. Polarized light images were combined into Stokes' 
parameters in AIPS, using standard formulae. 

Considerable effort was expended to estimate uncer- 
tainties. We accounted for several error sources, includ- 
ing zero-point and slope errors in the SYNPHOT calibra- 
tion, flat-fielding errors, and Poisson errors in the data 
and in modelling and subtraction of the galaxy. Errors 
were propagated by adding in quadrature. The cumula- 
tive error in the /?5T photometry is typically 3-4%, and 
is dominated by the flat-flelding and zero-point errors. 

Lightcurves were extracted from our data at three 
UV/optical wavelengths - 220 nm, 475 nm and 814 nm, 
and three X-ray bands (paper I) - 0.2-0.75 keV ("soft"; 
nominal 0.5 keV), 0.75-2 keV ("medium"; nominal 1.4 
keV), and 2-6 keV ("hard"; nominal 4 keV). The 2003 
March 31 and 2003 May 10 (part of a snapshot program 
led by D. Maoz) images were taken at a slightly longer 
wavelength; to place them on a common 220 nm flux scale 
we assumed ao — 0.6 (Perlman et al. 2001a). Finally, we 
have also made use of flux points from 2001-2002 HST 
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Table 1. Log of HST Observations 



Date 



Instrument & Bands 



Program 



2002-11-30 ACS F220W, F475W, F814W 9705 

2002-12-07 ACS F606W + POLVIS 9705 

2002- 12-22 ACS F220W, F475W, F814W 9705 

2003- 02-02 ACS F220W, F475W, F814W 9705 
2003-03-06 ACS F220W, F475W, F814W 9705 
2003-03-31 ACS F250W, F330W 9454 
2003-04-17 ACS F220W, F475W, F814W 9705 
2003-05-10 ACS F250W, F330W 9454 
2003-06-07 STIS F25QTZ 9474 



data which are discussed in more detail in Paper III. 

3. RESULTS 

Figure 1 shows the 2003 April 17 F220W image of the 
inner 4.5" of the M87 jet, as well as optical and X-ray 
lightcurves for the nucleus and HST-1 respectively. As 
can be seen, both components vary on timescales of a 
few months in the optical, with typical month-to-month 
variations being about 10-15% at 220 nm. Other regions 
do not vary significantly from month to month and so 
are not shown in Figure 1. For the nucleus, flux was 
extracted from 0.2" and 0.6" square boxes, which showed 
identical variability. In Figure 1 wc use the 0.2" box. 
For knot HST-l, we extracted flux from a 0.6" box. The 
extraction regions for our Chandra data were described in 
paper I. Our data conflrm Tsvetanov et al.'s (1998) claim 
of month-timescale nuclear optical variability. However, 
this is the flrst such report for any jot component. 

At 220 nm, > 95% of the flux in both varying regions is 
contained within a distribution consistent with the PRF. 
Thus both arc unresolved, and so < 0.02" (1.5 pc) in size. 
This implies light-travel times within a factor 3-20 (de- 
pending on the Dopplcr factor S = [T{1 — /3cos9)]~^) 
of the variability timescales (§4, Paper I). The nucleus's 
variable region is < 0.01" from the galactic center, with 
HST-l's varying region being 0.82" (64 pc) distant, loca- 
tions identical to ± 0.01" of those measured with Chan- 
dra (Paper I). We detect no motion in HST-l's flaring 
region; therefore, the flare likely occurred in the nearly 
stationary component at its upstream end (Biretta et 
al. 1999), rather than one that propagates at the veloc- 
ities seen in moving components within HST-1, namely 
Vapp = 6c (0.013"/ycar, easily detectable in these data). 

Figure 2 shows the F606W polarized light image. As 
can be seen, the nuc;leus is essentially unpolarized. HST- 
1 shows fractional polarization P = 0.46 at its flux peak, 
a maximum P = 0.68 (nearly the theoretical maximum 
for synchrotron radiation from an ordered magnetic field) 
at its upstream end ^ 0.72" from the nucleus, and a 
minimum P = 0.23 0.92" from the nucleus. The mag- 
netic field vectors in HST-1 are perpendicular to the jet 
direction, consistent with a shock. The alternative in- 
terpretation of a tightly wound helix would require cy- 
cles separated by < 0.02" to be consistent with the ob- 
served morphology. A much higher bulk F would also 
be needed (because of the longer path) to be consistent 
with Vapp = 6c (Biretta et al. 1999). The polarization at 
HST-l's flux peak is much higher than in 1995 (P = 0.14, 
Perlman et al. 1999). This is in line with the properties 
of BL Lac objects, which often have higher optical polar- 
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Fig. 1.— At top, we show the 2003 April 17 HST/ ACS F220W 
image of the inner 4.5" of the M87 jet. The image has been rotated 
so that the jet appears along the li-axis. The middle and bottom 
panels show lightcurves for the nucleus and HST-1. Some of the 
optical/UV lightcurves are multiplied by arbitrary factors to sep- 
arate them from the X-ray lightcurves. Flare events discussed in 
the text are enumerated. All error bars are at the la level. 



izations in high states (e.g., Hagen-Thorn et al. 1998). 

4. CONSTRAINTS FROM THE LIGHTCURVES 

To constrain the physics, we compute doubling and 
halving timescales T2, Ti/2, given in Table 2 for timeperi- 
ods referenced in Figure 1. We can compare these with 
predictions from models with E"^ losses, which are of 
the form t cc . If the emission is synchrotron ra- 
diation and synchrotron cooling dominates, rj = —0.5, 
(Kirk, Mastichiadis & Rioger 1998: Paper I). As shown 
in Bottcher et al. (2003), if synchrotron self-Compton 
(SSC) losses dominate, the dominant factor in d-y/dt is 
the synchrotron radiation energy density Usy, so that 
rj = {q — 4)/2 = -1 to -0.75 for electron spectral index 
q = 2.5 to 2 [^^(7) cx 7"*, which is related to a via 
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Fig. 2. — An image of the innermost 1.5" of the M87 jet, with 
polcirization vectors (magnetic field) superposed. The magnitude 
of the vector is proportional to the percent polarization, with a 
0.02" vector representing 100% polarization. Contours are shown 
at ( 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 4, 6, 10, 20) ADU/sec 



a={q-l)/2]. 

Table 2 shows several differences between the variabil- 
ity behavior of the nucleus and HST-1. The nucleus does 
not show strong energy-dependence in cither the X-rays 
{r|x,r^se = "0.24 ± 0.21 and rj^jaii = -0.07 ± 0.10) or 
optical {r]o,rise = —0.09 ± 0.43 and ijojaii = 0.52 ± 0.32, 
but the departure from zero in the latter measurement is 
due entirely to the lower significance of the fall in F220W 
during 2003 Feb. -Apr.; further, note that the drop accel- 
erated after monitoring in F475W and F814W ceased). 
By comparison, HST-1 does show energy-dependent be- 
havior in the X-rays, with r}x,rise = —0.62 ± 0.13 and 
tlx, fall = —0.49 ± 0.21. But in the optical, our data 
are less constraining, giving rjo rise = —0.35 ± 0.22 and 
riojaii = -0.51 ± 0.78. 

HST-l's behavior in the X-rays agrees with the predic- 
tions of a simple synchrotron model. But such a model 
cannot easily account for the comparable optical and X- 
ray variability timcscalcs. Two other models can explain 
such behavior. The first is that the flare was caused 
by adiabatic compression, followed by expansion on a 
dynamical timescale of Tdyn ~ 200 days. This would 
allow radiative losses dominate at high energies, where 
Tsyn < Tdyn, whilc cxpansion losses would dominate at 
low energies. The X-ray variability behavior then re- 
quires B K.2 mG. However, this model has difficulty 
explaining X-ray and optical increases and decreases that 
do not exactly coincide in time (Figure 1). Alternately, 
the X-ray fiarc could be triggered by shock compression 
(as suggested by the polarimetry, §3), with the optical 
emission representing the shocked plasma's downstream 
propagation. This model requires a ~ lOx stronger mag- 
netic field downstream of the shock to explain the similar 
0.5 keV and 220 nm variability timescales. 

The nearly energy-independent variability behavior of 
the; nucleus disagrees with simple synchrotron models, 
although some cannot be excluded formally because of 
the large error bars. One could account for the energy 
independence either by adiabatic compression and ex- 
pansion, or a helical trajectory, where flux would change 
with viewing angle (suggested by Urry et al. 1993 to ex- 



plain the nearly energy-independent variability of PKS 
2155—304 in 1991). But perhaps the most attractive 
model is that the varying region has R/c « 70 5 light- 
days, so that the light-crossing time controls the variabil- 
ity behavior (e.g., Chiaberge & Ghisellini 1999). In this 
case we would require _B ^ 20 (5^^ mG. 

5. COMPARISON TO BLAZAR FLARES 

The flares of the nucleus and HST-1 (Figure 2) resem- 
ble those of blazars, albeit on longer timescales. Can we 
explain blazar flares (t ~ 0.2—5 days; Ulrich, Maraschi & 
Urry 1997; Plan 2002), as highly beamed versions of these 
events? Timcscalcs Lorcntz-transform as r' oc t5~^ , so 
T = 5 days requires Smst/ ^blazar ~ 1/12 (also consistent 
with the luminosities seen in blazar flares, e.g., Giebels et 
al. 2002), while r = 1 day requires SMsr/^uazar ~ 1/60. 
For reasonable F, this requires Sm87 ~ 2, at the lower 
end of the range considered for HST-1 in Paper 1, and 
implies Suazar = 20 — 100. Consistency with HST-l's 
Vapp = 6c (Biretta et al. 1999), then requires F > 10 and 
^ = 15 — 19°, much tighter constraints than in previous 
works. A timescale r = 5 days then requires F = 10 
and = 1° if 5ms7 = 1-6 (^ Fm87 = 10, 0m87 = 19°), 
while T = 1 day requires F > 25 and 9 < 1°, plus smaller 
(^ii/ST- Unifled schemes usually constrain F and 6 less 
tightly (e.g., F = 5 - 30 and 61 < 10 - 30°, Urry & 
Padovani 1995). However, one can imagine a range of 
S or stronger magnetic flelds (e.g., Li & Kusunose 2000, 
Kataoka et al. 2000, Bottcher et al. 2003). 

Blazars can show either energy-dependent or energy- 
independent behavior, even within the same campaign. 
Often, energy-dependent behavior includes one band 
leading the other. With our data, any discussion of de- 
lays is premature, although there are possible indications 
in the data. Monitoring of M87 continues, and we will 
discuss this subject in a later paper. 

6. SPECTRAL VARIABILITY 

Do we also see spectral variability in these compo- 
nents? Figure 3 shows broadband spectra in three 

epochs, representing both high and low states. 

The nucleus has cto > 1 and aox > 1, in agreement 
with the flt value of Vbreak ~ 3 x 10^^ Hz found using 
historical data (Perlman et al. 2001b). But in all epochs, 
C(x < Cfox < cto, where ao,aox and are, respectively, 
the spectral indices in the optical, optical- X-ray and X- 
ray bands. Thus for standard synchrotron models, the X- 
ray emission must come from a component distinct from 
that responsible for the optical emission, since a single 
electron population cannot produce spectral hardenings 
blueward of Vbreak- The nucleus displays only modest 
spectral variability, with the 2002 December points (a 
low optical state) having flatter ao but steeper aox- This 
is consistent with a significant SSC contribution. 

By contrast, HST-1 has ao < aox < ctx in all epochs, 
with ao = (0.41,0.47,0.59) for epochs (2003 April, 2002 
December, 2003 February), similar to that seen in 1998, 
when it was 10 x fainter in optical (Perlman et al. 2001a). 
These values of ao are consistent with shock injection. 
The optical to X-ray spectrum steepens in low states and 
flattens in high states, as expected for synchrotron flar- 
ing. Interestingly, the 2003 April high state has a clearly 
harder spectrum than that seen in 2002 December. 
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Table 2. Flare Timescales 



ComporiGnt 


Band 


Dates 


lype (f lare ^) 




Doubling or 
4T5 nm 


Halving Timescale in Band (days) 
220 nm 0.5 keV 1.4 keV 


4 keV 


Nucleus 


X-rays 


June-July 2002 


Rise (2) 








135 ± 39 


608 ± 474 


83 ±20 


Nucleus 


X-rays 


March-May 2002 


Fall (1) 








78 ±11 


71 ±6 


67 ±9 


Nucleus 


Optical 


Dec. 2002-Feb. 2003 


Rise (2) 


317 ± 108 


206 ± 47 


256 ± 82 








Nucleus 


Optical 


Mar.-Apr. 2003 


Fall (2) 


81 ± 12 


83 ± 12 


274 ± 106 








HST-1 


X-rays 


Feb.-Apr. 2003 


Rise (2) 








223 ± 49 


118 ± 13 


61 ±9 


HST-1 


X-rays 


Nov.-Dec. 2002 


Fall (1) 








149 ± 55 


91 ± 17 


54± 10 


HST-1 


Optical 


Feb.-Apr. 2003 


Rise (2) 


359 ± 82 


414 ± 108 


243 ± 36 








HST-1 


Optical 


Dec. 2002 - Feb. 2003 


Fall (1) 


260 ± 133 


217 ±91 


139 ± 37 









Nucleus 



000.0 


























a 


















& 


A 




UU J Mp III ^ 










100,0 










2002 December 


















A 


2003 February (/2) 








I 


10.0 


















: 


1 .0 












B 






















a 


















A 




H 




0.1 








1 




1 


A 


* 




14 






1 b 


16 


17 




18 














log Frequency (Hz) 




















H=--T-1 










000.0 
























a 


a 






































a 


2003 April (x2) 










100,0 






A 




2002 December 




















2005 February (/2) 










10.0 




















1-0 












a 






















a 


















A. 




a 




0.1 














A 


4 




14 








16 


17 




15 





log Frequency (Hz) 

Fig. 3. — Spectral Energy Distributions for the nucleus (top) and 
knot HST-1 (bottom). Three epochs are shown; 2002 December, 
2003 February, and 2003 April. 
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7. FINAL THOUGHTS 

Here and in paper I we have found strong optical and 
X-ray variability in the nucleus of M87, and knot HST- 
1 in its jet. Is it possible to link the variability of the 
nucleus and HST-1 in any way? Perhaps the flaring in 
HST-1 was triggered by a density enhancement travelling 
down the beam of the jet. If we assume a near-constant 
apparent speed of Vapp = 6c for this material, its ejection 
would have occurred « 30—35 years ago. Such an episode 
could plausibly have caused a radio flare in the nucleus. 
Interestingly, DeYoung (1971) detected ^ 30% variations 
in the nuclear radio flux in 1969-1971. It is possible (but 
speculative) to link the flaring in the nucleus in ~ 1970 
to the current behavior of HST-1. 
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